Non Perturbative Renormalization Group study of reaction-diffusion processes 
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We generalize nonperturbative renormalization group methods to nonequilibrium critical phenom- 
ena. Within this formalism, reaction-diffusion processes are described by a scale dependent effective 
action, which flow is derived. We investigate branching and annihilating random walks with odd 
number of offsprings. Along with recovering their universal physics (described by the directed per- 
colation universality class) , we determine their phase diagrams and predict that a transition occurs 
even in three dimensions, contrarily to what perturbation theory suggests. 
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The description of many physical, chemical, and bio- 
logical systems resorts to nonequilibrium many-particle 
models. Their dynamical evolution is governed by a 
Markov process ruled by microscopic transition prob- 
abilities which violate detailed balance. The ensueing 
macroscopic behaviors are hence intrinsically much richer 
than in equilibrium. Many models exhibit phase transi- 
tions between stationary states, which have unveiled the 
existence of universal properties. The characterisation 
of the nonequilibrium universality classes constitutes a 
prevalent issue of statistical mechanics. In this respect, 
reaction-diffusion processes provide a simple and thus 
valuable framework to investigate nonequilibrium criti- 
cal phenomena. They generically describe diffusing par- 
ticles that can undergo different microscopic reactions 
P, Q • The most common universality class encountered 
is that of directed percolation (DP) It can be sim- 

ply achieved through a reaction-diffusion model endowed 
with the processes: 

2 A — > pair annihilation at rate A (1) 
A — > 2A generation of an offspring at rate a (2) 
A — > spontaneous decay at rate p (3) 

The competition between birth and death gives rise to 
a steady state where the density saturates at a constant 
value p. The system exhibits a continuous phase transi- 
tion between an active (fluctuating) state where p > 0, 
and an inactive (absorbing) state where the density van- 
ishes. In fact, a large variety of models, like the contact 
process , Reggeon Field Theory in particle physics , 
the dimer poisoning Q, the autocatalytic reaction mod- 
els 0, show the same critical behavior. 

Despite considerable successes, the general under- 
standing of nonequilibrium phenomena suffers from the 
lack of a systematic theoretical framework. Indeed, out 
of equilibrium, the absence of a free-energy functional 
from which to draw the probability distribution of steady 
states prevents from the use of well-established equilib- 
rium methods. Instead, the master equation should be 



solved, which is not feasible in general. Fortunately, it 
can often be cast into a path integral representation |9(, 
rooting a field theory that can be probed with a dynam- 
ical renormalization group (RG). The usual perturbative 
techniques consist in expansions in series of a parame- 
ter. They have allowed to achieve numerous advances in 
our theoretical understanding. Nevertheless, as they rely 
on the smallness of a parameter, their validity remains 
confined to the vicinity of critical dimensions or to weak 
coupling regimes. 

We propose in this Letter a powerful alternative, based 
on a nonperturbative RG procedure (NPRG) |l0lllllll2| . 
This method avoids the previous pitfalls as it does not 
intrinsically depend on the magnitude of a parameter. 
Moreover, the microscopic interactions can be traced all 
along the NPRG flow, and related to the long-range 
properties, thus giving access to non universal quanti- 
ties, such as equations of state and critical temperatures 
(lij. In this work, we generalize the NPRG methods to 
nonequilibrium systems and derive very generic, model- 
independent RG flow equations for reaction-diffusion pro- 
cesses. We exploit these equations to investigate branch- 
ing and annihilating random walks (BARW) with odd 
number of offsprings. We first recover their universal 
physics which belong to the DP universality class. Then, 
linking the density p of stationary states with the micro- 
scopic rates enables us to push further a recent analysis 
by Cardy and Tauber Q . We predict in particular that 
a phase transition can occur in three dimensions, con- 
trarily to what was conjectured. 

The NPRG formalism — the effective average action 
method |ld| — is a continuous implementation, on suc- 
cessive momentum shells, of Wilson's block-spin concept. 
For equilibrium systems, it consists in building a se- 
quence of effective actions T^, at the running momentum 
scale k, interpolating between the microscopic Hamilto- 
nian H of the system and its Gibbs free energy T. This 
is achieved by adding to the partition function Z a mass- 
like term AHk that suppresses the propagation of the 
low-energy modes (with wave number q < k) while not 
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affecting the high-energy ones (q > k). The correspond- 
ing free energy \ogZ k , or its Legendre transform 1^, thus 
only encompasses the large momentum fluctuations. At 
the scale k = A, A -1 corresponding to the underlying 
lattice spacing, no fluctuation has yet been taken into 
account and Tk=A = H, while at k = 0, all fluctuations 
have been integrated out and Tk=o = T. An exact RG 
equation for the flow of T\ is easily derived and underlies 
the NPRG formalism [lOj. We now generalize this for- 
malism to out of equilibrium systems, for which partition 
function and Gibbs free energy no longer exist. Despite 
this lack, a path integral formulation of the dynamics 
can be set up, commonly involving an auxiliary response 
field (j>, along with the usual field tj>. It allows to define a 
generating functional Z for the correlation functions, and 
one, F, for the one-particle irreducible correlation func- 
tions. These functionals constitute the building-blocks 
of an out of equilibrium NPRG. We introduce in Z a 
scale-dependent term, analogous to AiJ^: 



Z k [J,J] = / D<j>Dj> e 



-S[0,0]-A5 fc [0,0]+/ J<p+f J<p 



(4) 



where 



d d p duj 
(27r) d 2^ 



1>(-p,-uj)Rk(p)<P T (p,Lu), (5) 



$ denoting the two-component vector (0, (f>) . The 2x2 
matrix Rk is symmetric with zeros on its diagonal, and 
a mass-like function Rk{p) off its diagonal. The selec- 
tion of modes is ensured by imposing that, on the one 
hand, Rk(p) ~ k 2 for p <C k, that is R k acts as an ef- 
fective mass — an IR cut-off — for slow modes. On the 
other hand, one imposes that Rk(p) — > for p ^> k, so 
that the rapid modes remain unaltered. In this work, we 
use Rk(p) = (k 2 — p 2 )6(k 2 — p 2 ) which gives sim- 
ple analytical expressions for the momentum integrals. 
We emphasize that no cutoff is necessary for time, since 
the corresponding integrals are regular and the frequency 
dependence can be integrated over (see below). The ef- 
fective action Tk is defined through the appropriate Leg- 
endre transform Tk + logZ k — J Jip + Jip — AS k (ip,ip), 
where ip (resp. ip) is the mean value of <p (resp. </>), de- 
fined as the derivative of Zk w.r.t. J (resp. J). The 
mass term ASk ensures that the proper limits Tk=A = S 
and Tk=o — T are recovered. The exact flow equation 
of Tk is then easily derived, similarly to the equilibrium 
calculation [Io| : 



dkTk(^,^) = 



• (2) 

where T k is the field dependent matrix of the second 
functional derivatives of T, and Tr stands for matrix trace 
over internal indices and integration over the internal mo- 
mentum and frequency. Eq. © is a functional equation 



±Tr{[£W + R k ]- x d k R k }, (6) 



which obviously cannot be solved exactly. To handle it, 
one has to truncate Tk- Since the critical physics corre- 
sponds to the long distance (q — > 0), long time (w — » 0) 
limit, a sensible truncation consists in expanding Tk in 
powers of gradients [Toj and time derivatives. Retaining 
only the leading order in derivatives, the functional Tk 
for reaction-diffusion systems reads: 

r*0M) = J d d rdt[Zk^(dt-DkW 2 )^ + Vk(^,^)}. (7) 

In principle, the three functions V k , Zk and Dk depend 
on the fields, and their functional form is dictated by the 
symmetries of the underlying dynamics. We study in the 
following the leading order approximation that consists in 
neglecting the field-dependence of Z k and Dk 01- Thus, 
we only consider their fc-dependence, which accounts for 
the critical anomalous scalings of the fields, and of time. 
Indeed, the anomalous dimension rj of the fields is de- 
fined such that at criticality, ip and i/5 scale as £;( d + r ')/ 2 . 
This definition induces that 77 = — dlnZk/d\nk. The ex- 
ponent (z — 2) embodies the anomalous scaling between 
time and space in the critical scaling regime jlfij . It fol- 
lows that 2 — z = —dlnDk/dlnk. We emphasize that the 
previous ansatz Q should qualitatively well encompass 
the physics of the model, and moreover provide accurate 
estimates of "static" quantities, related to the potential 
part. However, this level of approximation is not ex- 
pected to give very accurate dynamical exponents, even 
less that they are large. Refining these exponents would 
require to include the field dependence of Zk and Dk, 
and as a further step, to enrich the derivative content of 
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We can now establish the generic flow equations for 
reaction-diffusion processes. The flow of Vk is drawn from 
Eq. iJBJ, evaluated at a uniform and stationary field con- 
figuration. After integrating over u, one gets: 



dkVk 



1 

2Zk 



d d p 
(2^p 



dkRk 



dlv k d 2 v k 

hk( P ) 2 



(8) 



where h k (p) = Z k D k p 2 + R k (p) + d 2 -^V k . The flow 
of Zk (resp. Dk) is obtained taking the q 2 (resp. v) 
part of the flow of <9| ^T k {p, v). These flows depend on 
the fields, but the resulting exponents should be identi- 
cal for any field values. However, this property breaks 
down when any truncation is performed |I3j. Since, in 
the reaction-diffusion models, <f> vanishes, and the av- 
erage of 4> reaches a uniform and stationary configura- 
tion, it is sensible to define the exponents at ip = and 
ip{p,v) = ipokS d (p)S(v). Derivating twice Eq. © w.r.t. 
i}){—q,—v) and ip{q,u) yields, after integrating over the 
internal frequency: 
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d k d^T k (q, v) 



2Z ^Tpip^ k ^Pilnpil> 



Ad 

a p 



dkRk{p) 



{2n) d [Z k w + h k {p) + h k (p-q)f 



v k y 



d> 
(27r) d 



dkRk{p) 



hk{p) + h k (p- q) 



h k (p)h k (p - q)[Z k \v + h k (p) + h k (p - q)Y x 



h k (p) 2 [Z k iv + h k (p) + h k {p- q)} 



(9) 



As a first step, we show that these equations allow to 
recover the universal physics of DP and more precisely 
to compute the three usually quoted independent critical 
exponents z, v and [3. v describes the scaling behavior 
of the correlation length near criticality, (3 that of the 
order parameter. A held theoretical action S[<f>, (b] can 
be derived from the microscopic processes {Q, @ and 
(J3J, upon the introduction of two fields, 4> and (f>. After 
their rescaling and a shift of 4>, one gets : 



d d rdt [<b(d t -DV 2 



l2 ) + X<b 2 $ 2 



(10) 



which is invariant under the simultaneous changes 
4>{x,t) — > —cf)(x,—t) and (b(x,t) — > —cj>(x,—t). The 
generic term of the effective potential V k ('tp, ip) allowed by 
this symmetry reads a a p(^> a ^ + (— l) a+l3 ^p a 'ifj 13 ), which 
can be expressed in terms of the two invariants x = ipip 
and y = ip—ip. Thus V k can be conveniently parametrized 
by U k (x,y) = Vk^^ip). By numerically integrating the 
flows and ©, we find that for a fine tuned initial 
(fx — a), corresponding to criticality, the (dimensionless) 
effective potential flows to a fixed function. The expo- 
nents are calculated in the vicinity of this fixed solu- 
tion. We first implement the well-known local poten- 
tial approximation (LPA), that is Z k = D k = 1. The 
static exponent v is very accurate in all dimensions, see 
Table whereas /3, deduced from the scaling relation 
[3 = l/2(d + r/), is more imprecise since t] = at this 
order. We then study the leading order (LO) in deriva- 
tives: Z k and D k now run with k. The critical expo- 
nents are all the more accurate that rj and (z — 2) remain 
small. They thus deteriorate as the dimension lowers, 
since rj and (z — 2) then become large. As usual with this 
method, more accurate determinations of rj and z rely 
on a proper description of the momentum and frequency 
dependence of the two-point correlation function, which 
would require to implement a richer ansatz for T k |l8j . 
as already emphasized. 

We finally come to the most important part of this 
work: the study of non universal features of the BARW 
defined by the set of the two processes and (J2J 
(// = 0). Mean field theory predicts the density to reach 
a finite saturation value for any positive birth rate a, and 



d 




(a) LPA 


(b) LO 


(c) MC 




V 


0.584 


0.548 


0.581(5) 


3 


P 


0.872 


0.782 


0.81(1) 




z 


2 


1.909 


1.90(1) 




V 


0.730 


0.623 


0.734(4) 


2 


13 


0.730 


0.597 


0.584(4) 




z 


2 


1.884 


1.76(3) 




V 


1.056 


0.888 


1.096854(4) 


1 


p 


0.528 


0.505 


0.276486(8) 




z 


2 


1.899 


1.580745(10) 



TABLE I: Critical exponents of DP. (a) and (b): NPRG calcu- 
lations from this work, at the local potential approximation, 
resp . at leading order, (c): from Monte Carlo simulations 



thus the system to be always in the active state. Never- 
theless, simulations |20j have brought out the existence 
of a phase transition in this model in d = 1 and d = 2, 
invalidating mean field conclusions. Cardy and Tauber 
have addressed this question through a field theoretical 
perturbative study |f l| and provided a partial answer rel- 
evant for small reaction rates X/D and a/D, and in the 
vicinity of d = 2. They have shown that fluctuations 
can indeed induce, for a non zero critical birth rate cr c , 
a dynamical transition to an absorbing state (belonging 
to the DP universality class). The underlying reason for 
the existence of the transition is that the combination 
of (JTJ and generates under renormalization the one- 
particle spontaneous decay J2J). This process allows the 
existence of an absorbing state if its renormalized rate 
/ifj(cr, A) renders /«r — <jr positive (see Eq. llTJ|) . In d = 2, 
this defines a narrow absorbing region below a critical 
transition curve given by: 



a c (X)/D oc e 



-AttD/X 



(11) 



For d > 2, perturbation theory breaks down. However, 
since the transition curve Eq. Ijlljl is already infinitely 
flat for small X/D in d = 2, and since fluctuations are 
expected to be suppressed when the dimension increases, 
it has been suggested in 01 thai no transition should 
exist in d = 3. 

We re-investigate this question using NPRG, which is 
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neither restricted to small X/D nor to a specific dimen- 
sion. We can thus explore the phase diagrams of BARW 
in d = 2 and d = 3. They are displayed in Fig. QJ. 
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FIG. 1: Phase diagrams of BARW in d = 2 and d = 3. 

In d — 2, we find two regimes. For small X/D (typically 
below 4) the transition curve is very accurately fitted 
by an exponential with a coefficient —11.8, in agreement 
with Eq. C3J at a 6% level. For large X/D, a c {X)/D 
grows linearly opening up a large absorbing phase. In 
d = 3, we also find two regimes, one with a linear critical 
curve, and one without transition, separated by a thresh- 
old value X/D ~ 30. Let us now give simple arguments 
highlighting why a transition must exist in d = 3, that 
perturbation theory misses. 

When D is large compared with A and a, the mech- 
anism of particle destruction is driven by the pair anni- 
hilation occuring when wandering particles meet. While 
in d = 2 all the particles can be destroyed through this 
mechanism, it is not efficient enough to do so in d = 3. 
We thus expect only an active phase at large D in d = 3. 
This corresponds to the result of the perturbation the- 
ory, which is indeed valid at small X/D. On the other 
hand, in the small D regime where perturbation theory 
fails, the mechanism of destruction is dominated by "self- 
destruction" , that is particles annihilate with their own 
offsprings before they can diffuse. This process no longer 
requires the random encounter of particles and thus does 
not depend on the dimension. We therefore expect a 
transition to occur at small D both in two and three 
dimensions. These arguments corroborate the phase di- 
agrams obtained from the NPRG analysis. Though, re- 
lying on the ansatz (0, we cannot provide a quantitative 
estimate of the accuracy of the critical rates in Fig. ^ 
the previous arguments reinforce the NPRG ability to 
capture the non universal physics of the model over the 
different diffusive regimes. Moreover, we have performed 
numerical simulations, which confirm the existence of a 
DP phase transition in d = 3 for large X/D, and hence 
support our prediction (extensive simulations to quanti- 
tatively check the phase diagrams arc in progress 21]). 



Note that in simulations |2(j, the implementation of 
reaction involves the creation of offsprings on neigh- 
boring sites, contrarily to the on-site birth process under- 
lying field theoretical approaches (both |l4j and ours). 
Though this should not alter the behavior of the model 
when D is large, it becomes crucial at small D since 
the "self-destruction" mechanism is then drastically sup- 
pressed. This is probably why no transition is found in 
d = 3 in these simulations. 

In the present work, we have implemented for the 
first time NPRG methods in nonequilibrium statisti- 
cal physics, and derived very generic flow equations for 
reaction-diffusion models. This method has allowed us 
to recover the universal behavior of directed percolation 
and moreover to compute the phase diagrams of BARW. 
We have shown that a phase transition exists in d = 3 
in these systems contrarily to what suggests perturba- 
tion theory. The formalism developed here offers an as 
efficient as promising mean of investigation of nonequilib- 
rium physics. This could give some valuable insight into 
the perennial Kardar-Parisi-Zhang equation, describing 
roughening transitions in growing interfaces. 

We wish to thank H. Chate for providing us codes and 
assistance for simulations and F. Van Wijland and U. 
C. Tauber for helpful discussions. The LPTHE is Unite 
Mixte du CNRS, UMR 7589. 
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